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Abstract 

This paper studies the risk-adjusted optimal timing to liquidate an option at the prevailing 
market price. In addition to maximizing the expected discounted return from option sale, 
we incorporate a path-dependent risk penalty based on shortfall or quadratic variation of the 
option price up to the liquidation time. We establish the conditions under which it is optimal 
to immediately liquidate or hold the option position through expiration. Furthermore, we study 
the variational inequality associated with the optimal stopping problem, and prove the existence 
and uniqueness of a strong solution. A series of analytical and numerical results are provided 
to illustrate the non-trivial optimal liquidation strategies under geometric Brownian motion 
(GBM) and exponential Ornstein-Uhlenbeck models. We examine the combined effects of price 
dynamics and risk penalty on the sell and delay regions for various options. In addition, we 
obtain an explicit closed-form solution for the liquidation of a stock with quadratic penalty 
under the GBM model. 
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1 Introduction 


For decades, options have been widely used as a tool for investment and risk management. As 
of 2012, the daily market notional for S&P 500 options is about US$90 billion and the average 
daily volume has grown rapidly from 119,808 in 2002 to 839,108 as of Jan 201 sB. Emp irical s tudie s 
on options returns often assume that the options are held to maturity (see Broadie et al. ( 20091 ) 
and references therein). For every liquidly traded option, there is an embedded timing flexibility 
to liquidate the position through the market prior to expiry. Hence, an important question for 
effective risk management is: when is the best time to sell an option? In this paper, we propose 
a risk-adjusted optimal stopping framework to address this problem for a variety of options under 
different underlying price dynamics. 

In addition to maximizing the expected discounted market value to be received from option sale, 
we incorporate a risk penalty that accounts for adverse price movements till the liquidation time. 
For every candidate strategy, we measure the associated risk by integrating over time the realized 
shortfall, or more generally its transformation in terms of a loss function, of the option position. As 
such, our integrated shortfall risk penalty is path dependent and introduces the trade-off between 
risk and return for every liquidation timing strategy. 

Under a general diffusion model for the underlying stock price, we formulate an optimal stopping 
problem that includes an integral penalization term. To this end, we define and apply the concept 
of optimal liquidation premium which represents the additional value from optimally waiting to 
sell, as opposed to immediate liquidation. As it turns out, it is optimal for the option holder to 
sell as soon as this premium vanishes. This observation leads to a number of useful mathemati¬ 
cal characterizations and financial interpretations of the optimal liquidation strategies for various 
positions. 

We first identify the conditions under which it is optimal to immediately liquidate or hold 
the option position through expiration. The investigation of the non-trivial liquidation strategies 
involves the analytical and numerical studies of the inhomog eneous variat ional inequality associ¬ 
ated with the optimal stopping problem. In a related work, Surya (2 0121 ) examines the solution 
structure for a finite maturity optimal stopping problem under Levy processes with a running 
cost and other features, and an inhomogeneous variational inequality also arises from the asso- 
ciated partial i ntegr o-differential free-boundary problem. In the context of asset management, 
Davanik and Egami ( 2012h study a perpetual optimal stopping problem with a running cash flow 
generated from dividend and coupon payments, and they solve a time-independent inhomogeneous 
variational inequality. For the variational inequa lities in our liqui dati on problems, we prove the 
existence and uniqueness of a strong solution a la Bensoussan and Lionsl ( 1978h (see Section [71 be¬ 
low) u nder general conditions applicable to both geomet ric Brownian motion (GBM) (se e I Merton 
( 19731 )) and exponential Ornstein-Uhlenbeck (OU) (see Ornstein and Uhlenbeck ( 193ol ) ) models 
for the underlying dynamics. We also provide some mathematical characterizations and numerical 
examples of the optimal liquidation strategies for stocks, calls, puts, and straddles. 

The incorporation of the risk penalty gives rise to optimal liquidation strategies that are dis¬ 
tinctly different from the unpenalized case. For instance, if the option’s Delta (derivative of the 
option price with respect to the underlying price) is of the same constant sign as the excess return 
of the underlying, then it is optimal to hold the option till maturity when there is no risk penalty 
(see Prop. m- This applies to the case of a call option (resp. a put option) if the investor is bullish 
(resp. bearish) on the stock. However, under risk penalization the investor may find it optimal to 
liquidate a call (resp. a put) early even in the bullish (resp. bearish) scenario (see e.g. Prop. 13.11) . 


x See http://www.cboe.com/micro/spx/introduction.aspx 
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Furthermore, the shape of the optimal liquidation region depends significantly on the risk penalty. 
We show that higher risk penalization coefficient always reduces the delay region, which intuitively 
means that the investor is more likely to sell earlier. Moreover, in some cases the optimal delay and 
sell regions can exhibit some interesting structures, such as disconnectedness (see Figures [3] and 
[6|) . These analytical results are significantly facilitated by the properties of the optimal liquidation 
premium (see Theorems 12.11 and 12 .411 . 

Our path-dependent risk penalization model can also be viewed as an alternative way to in¬ 
corporate the investor’s risk sensitivity in option liquidation/exer cise timing problems, a s com- 
pared to the utility maximization/indifference pricing approach ( Henderson and Hobson) . 1201 ll : 
Leung and Ludkovski . 2012 : Leung et al. . 20121 ). On the other hand, Leung and Ludkovskil (1201 il l 
investigate the optimal timing to buy equity European and American options without risk penalty 
under incomplete markets, w here the invest or is assumed to select risk-neutral pricing measure 
different from the market’s. Leung and Liu 1 20131 ) also discuss the timing to sell an option under 
the GBM model without any risk penalty, which is a special example of our model. 

As is well known, the concept of risk measures based on s hortfall risk has been applied t o 
many portfolio optimizati o n problems; see Artzner et al. ( 1999 1: Rockafellar and Urvasev ( 2000l ) : 


Follmer and Schiedl ( 2002 ): Follmer and Schied ( 20041 ) : lilhg 


an et al. ( 120051 ). and references therein. 


Our model applies this idea to options trading as a path-penalty associated with each liquidation 
strategy. As a variation of the shortfall we also introduce a risk penalty based on the quadratic 
variation of option price process. In particular, we obtain an explicit closed-form solution for the 
liquidation of a stock with quadratic penalty under the GBM model (see Theorem 15.11) . Through 
examining the optimal liquidation premium, we also compare the liquidation s trate gies for calls 
and puts under the shortfall-based and quadratic risk penalties. Fo rsyth et ajJ (120121 ) also adopt 
the mean-quadratic-variation as a criterion for determining the optimal stock trading strategy in 
the presence of price im pact. _ 

The recent paper bv lMacLean et al. (2013() considers a discrete-time portfolio optimization prob¬ 
lem with a convex loss function that accounts for the shortfall of the wealth trajectory from a bench¬ 
mark. While we consider the problem of optimal liquidation of stocks and optio ns, their investiga¬ 
tion focuses on the optimal capital growth or Kelly strategy. On the other hand, iFrei and West ray 
( 201bI ) study the optimal liquidation of a stock position subject to temporary price impact. Specifi¬ 
cally, they minimize the mean and variance of the order slippage with respect to the VWAP (volume 
weighted average price) as the benchmark. These papers adopt the stochastic control approach to 
solve for the optimal position over time, whereas our problems concern only the optimal timing to 
liquidate. 

The rest of the paper is organized as follows. In Section [2j we formulate the optimal liquidation 
problem for a generic European claim in a diffusion market. In subsequent sections, we focus on 
the liquidation of a stock or an option under the GBM and exponential OU models. In Sections 
[3] and [I] we study the optimal liquidation timing with a shortfall risk penalty. In Section [5j we 
conduct our analysis with a quadratic variation risk penalty. Section [6] concludes the paper. In 
Section [71 we discuss the existence of a strong solution to the variational inequality as well as the 
probabilistic representation satisfied by the optimal liquidation premium. 


2 Problem Overview 

In the background, we fix a probability space (11, T , P), where P is the historical probability measure. 
The market consists of a risky asset S and a money market account with a constant positive interest 
rate r. The risky asset price is modeled by a positive diffusion process following the stochastic 
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differential equation 


dS t = n{t, S t )S t dt + a(t, S t )S t dW t , S 0 = s, (2.1) 

where W is a standard Brownian motion under measure P and s > 0. Here, the deterministic 
coeffi cients u(t , s) and a (t.s) are assumed to satisfy standard Lipschitz and growth conditions 
( Karatzas and Shreve. 1 9911 . §5.2) to ensure a unique strong solution to (12.11) . We let F = [fFt)t >o 
be the filtration generated by the Brownian motion W. 

Let us consider a market-traded European option with payoff Ii(St) on expiration date T written 
on the underlying asset S. If the Sharpe ratio A (t,s) := r satisfies the Novikov condition: 

E{exp(J 0 T A A 2 (u, S u ) du )} < oo, the density process 


dP 


= exp 


T t 


ft 

A 2 (u,S u )du + / A (u,S u )dW v 
Jo 


0 <t<T, 


( 2 . 2 ) 


is a (P,F)—martingale (see iKaratzas and Shreve (|199lh . Prop. 3.5.12). This defines a unique 
equivalent martingale (risk-neutral) measure Q, and the market price of the option is given by 


V(t,s) = E M {e-^-^HST)} , ( t,s ) € [0,T] x 


(2.3) 


The shorthand notation E^, s { •} = E{-|5* = s} denotes the conditional expectation under Q. Note 
that the market price function V(t,s) does not depend on the drift function fi(t,s). 

Observing the stock and option price movements over time, the investor has the timing flexibility 
to sell the option before expiry. While seeking to maximize the expected discounted market value of 
the option, we incorporate a risk penalty that accounts for the downside risk up to the liquidation 
time. Specifically, we define the shortfall at time t by 


lit, S t ) = (■ m-V(t,S t ))+, 


(2.4) 


where m > 0 is a constant benchmark set by the investor. Then, the risk penalty is modeled as a 
loss function of the shortfall, denoted by St)). Here, the loss function ip : M + —>■ M is assu med 

to be increasing, convex, continuously differentiable, with ip( 0) = 0 (see e.g. ( Follmer and Schiedl . 
2004, Chap 4.9)). As a result, the investor faces the penalized optimal stopping problem 


J a (t,s) = 


sup 
-€7) T 


E 


't,S 


,—r(r— 


4 V(r, S T )-a J e~ r ( u ~Vip ((m - V(t, S t )) + ) du 


(2.5) 


where a > 0 is a penalization coefficient and % t is the set of E-stopping times taking values in 
[f,T], 

Unless otherwise noted, our analysis applies to a general loss function ip satisfying the conditions 
above. Here, let us give an example to visualize the penalization mechanism. For instance, one can 
set the benchmark to be the initial option price, and take ip(£) = l. Then, the penalty term amounts 
to accumulating the (discounted) area when the option is below its initial cost. We illustrate this in 
Figure[l] Notice that the realized shortfall stays flat when the option price is above the benchmark, 
and continues to increase as long as the option is under water. Other viable specifications include 
the power penalty ip(£) = £ p , p > 1, and the exponential penalty ip(£) = exp(y^) — 1, 7 > 0 , and 
more. 

In order to quantify the value of optimal waiting, we define the optimal liquidation premium by 
the difference between the value function J a and the current market price of the option, namely, 


L a (t,s) := J a (t,s) -V(t,s). 


( 2 . 6 ) 
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Figure 1: The realized shortfall (dashed) based on a simulated price path (solid) of a European call option 
under the GBM model, with parameters Sq = 100, r = 0.03, = —0.05 and a = 0.3, K = 100, T = 1, 
a = 1. The benchmark m is the initial call option price. 


Alternatively, the optimal liquidation premium L a can be interpreted as the risk-adjusted expected 
return from a simple buy-now-sell-later strategy. 

Denote the discounted penalized liquidation value process by 

Y u = e~ ru V(u, S u ) -a [ U e~ rt ^{(m - V{t, S t )) + )dt. 

Jo 


In order to guarantee the existence of an optimal stopping time to problem (12.51) . we require that 
E{sup 0 < u < T Y u } < oo. For a European call option, the option value V(t,St ) is dominated by the 
stock price St, while the put option price is bounded by the strike price. Consequently, for any 
linear combination of calls and puts, it suffices to impose E{sup 0 < u <r S u } < oo. We also require 
that P{min 0<t< £ S) > 0} = 1, which means that the asset price stays strictly positive before any 
finite time t a.s. Then by standard optimal stopping theory (j Karatzas and Shreve . 1998, Theorem 
D.12), the optimal liquidation time, associated with L(t,s), is given by 


t* = inf{ u € [t, T} : L a (u, S u ) = 0 }. 


(2.7) 


In other words, it is optimal for the investor to sell the option as soon as the optimal liquidation 
premium L a vanishes, meaning that the timing flexibility has no value. Accordingly, the investor’s 
optimal liquidation strategy can be described by the sell region S and delay region T>, namely, 

5 = {(f, s) € [0, T] x M+ : L a (t, s) = 0}, (2.8) 

V = {(i, s) € [0, T] x M+ : L a (t, s) > 0}. (2.9) 


Our framework can be readily applied to the reverse problem of optimally timing to buy an 
option. This amounts to changing the sup to inf in L a . In this paper, we shall focus on the 
liquidation problem. 
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2.1 Analysis of the Optimal Liquidation Premium 

Theorem 2.1. Given the underlying price dynamics in m, the optimal liquidation premium 
admits the probabilistic representation 

L a {t,s)= sup Et iS ( f T e- r ^G a (u,S u )du\ , (2.10) 

tG71, t Ut J 

where we denote 

G a (t,s ) := (n(t,s) - r)sV s (t,s ) - aif ((m - V(t, s)) + ) . (2.11) 

Proof. Applying Ito’s formula to the market price in (12.311 , we get 

{e~ r ^V(T, 5 r )} - V(t, s ) = E t>s { jT e“ r S u ) - r)S u V s (u, S„)duj . 

Substituting this into the optimal liquidation premium in (|2.6D gives 

L a (t,s) = sup E tjS | f e - r (“-0 S u ) - r)S u V s (u,S u ) - aip((m - V(u, S„)) + )] du\. 

re7i, T l Jt ^ - v - s J 

= G a (u,Su ) ill ED 


□ 

We shall call G a (t,s ) in (12.111) the drive function. We observe that it depends on the Delta 
V s = ^ of the option and the penalty coefficient a reduces the drive function for every (f, s). Many 
properties of the optimal liquidation premium L a can be deduced by studying the drive function. 


Proposition 2.2. Let t € [0, T] be the current time. If the drive function G a (u,s) is positive, 
V(it, s) € [t, T ] x M + , then it is optimal to sell at maturity, namely, t* = T. If the drive function 
G a (u,s ) is negative, M(u,s ) € [t,T] x R + , then it is optimal to sell immediately, namely, r* = t. 

Proof. We observe from the integral in (12.101) that if the drive function G a is positive (resp. neg¬ 
ative), V(u,s) G [t,T} x E + , then we can maximize the expectation by selecting the largest (resp. 
smallest) stopping time, namely, t* = T (resp. r* = t). □ 

In particular, if V s (t,s ) and — r) are of different signs V (t,s), then the drive function 

G a is always negative, so it is optimal to sell immediately. Proposition 12.21 can also be applied to 
the perpetual case if we set T = oo. In general, the delay region always contains the region where 
the drive function is positive, namely, 


{G a > 0} C {L a > 0}, 


( 2 . 12 ) 


see e.g. (j Oksenda l and Suleml 120051 . Prop. 2.3). Intuitively, this means that if G(t,s ) > 0, then 


the investor should not sell immediately since an incremental positive infinitesimal premium can 
be obtained by waiting for an infinitesimally small amount of time. 

In addition, we can infer from ()2.10j) the ordering of optimal liquidation premium based on the 
drive function. 
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Corollary 2.3. Consider two options A and B, along with two penalty coefficients a a and as 
respectively. If the drive function of A dominates that of B, i.e. G°f^(t,s) > G°^(t, s),V(f, s) € 
\0,T] xl + , then the optimal liquidation premium for A, L a , A , dominates that for B, Lf B (t,s), i.e. 
L“ A (f,s ) >L% B (t,s),\/{t,s) € [0, T] x R+. 

The corollary allows us to compare the liquidation timing of different penalties. For example, 
for 0 < < « 2 ) we have G ai (t,s ) > G a2 (t,s ) for the same option. It follows from (12.71) and 

Corollary 12.31 that the optimal liquidation time with penalty aq is later than that with penalty « 2 - 
In general, a variety of delay and sell regions can occur depending on the underlying dynamics 
and option payoff. Next, we give sufficient conditions so that the delay region is bounded. 

Theorem 2.4. Let T < oo and S be time homogeneous. Then, the delay region is bounded provided 
that 

(i) 3c> 0 s.t. G a (t , s) < c for every (t, s) €E [0,T] x R + ; and 

(ii) there exist constants b, k > 0 such that G a (t,s ) < —b in [0,T] x [k, oo). 

Proof. Step 1. We find a function L(t,s ) that dominates L a (t,s ) and is decreasing in both t and 
s. To this end, we define 

G a (s) := ma x{G a {t,f) : (t,f) € [0,T] x [s,oo)}, 

L{t,s):= sup E m ( [ e-^G« {Su )du, 

T€Tt,T Ut 

By construction G a : [0, T] x R + —x R is constant in t and decreasing in s. It also satisfies conditions 
(i) and (ii). Consequently, using the time homogeneity of S, we have, for t > t', 

L(t, s) = sup E 0jS { [ e~ ru G a (S u ) du\ < sup E 0 , s ( [ e~ ru G a {S u ) du\ = L{t ', s). 
re7o,T-t (jo J r&r 0T _ t , l JO J 

Hence, L(t,s ) is decreasing in t. Moreover, since G a is decreasing in s, we have, for s' > s, 

L(t, s') = sup E ty j [ e - r ( u -^G a (S u )du 

T£Tt,T (ii 

= sup E^j, { f e~ r ( u ~^G a (S u + s' — s) du 

T&T t ,T ’ Ut 

< sup E t)S ( f T e- r ^d Q (S u )du\ =L(t,s). 

T£Tt,T Ut J 

Therefore, L[t,s ) is also decreasing in s. 

Since by definition G a dominates G a , Corollary 12.31 implies that L a (t, s ) has a bounded support 
as long as L(t, s) has a bounded support. Henceforth, we can assume without loss of generality 
that L a is decreasing in both variables t and s and that G a is time homogeneous and decreasing 
in s. In particular, we denote G a (s) = G a (t,s). 

Step 2. We prove that for every t > 0, there exists s < oo such that L a (t,s) = 0 for every 
s > s. Since L a (t,s ) is decreasing, it is equivalent to show that there exists no f G (0, T] s.t. 
L a (t, s) > 0 for 0 < t < t and s € M + . To this end, let’s suppose that such a time t exists. In other 
words, t* = inf{t < u <T : L(u , S u ) = 0} > t. 
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Now, we show that this leads to a contradiction. Fix t € [0 ,t). Condition (ii) means that there 
exists k s.t. G a (s ) < — b < 0 in [fc, oo). For s > k, we let Tfc := inf {it > t : S u < k}. Since S has 
continuous paths, we have Tk > t. Define 

{ f r * Arfc 

b J e ~ r (. u ~ t ) du 

= - b (l - E m |e- r ( r * ATfc -*)}) , (2.13) 




where c is the upper bound of G a in condition (i). Next, taking s t oo yields that P t,sijk < T) {. 0, 
while {e _r ( T * Ar,I_t )} < since r* > t > t a.s. Therefore, we obtain 




cPt, s (r fc < t *) 
r( 1 -E t)S { e - r ( r * Ar fc- t )}) 


-> 0 . 


As a result, for a sufficiently large s > k, we get b > /3(t,s), which implies that JC(t,s) < 0 (see 

TO). 

Next we consider the difference 


L“(t, S ) - /C(t, s) < E m { [ T e~ r ^G a (S u ) du - -e- r(r '=- t) l {Tfc < r * } 

( JT*/\Tfr V 

< e r, E t , {he- rr ' AT ‘ - e~ rT *) - ^-^l {T .< n , ) } 

=-^E,„ < 0. 

This means that L a (t, s ) < /C(t, s) < 0. This contradicts the assumption L a (t, s ) > 0. 

Step 3. It remains to show at time 0 that 3 s > 0 such that L“(0, s) = 0 for every s > s. Let 
t € [0, T] and consider for every t € [0, T + f] the optimal stopping problem 

L a (t,s):= sup E tjS ( f e~ r{u - t) G a {S u )du 



The time homogeneity of S yields that L a (t,s ) = L a (0,s). Now we apply Step 2 and conclude 
that there exists s > 0 such that L (t, s) = 0 for every s > s. Hence, the delay region is bounded 
above. □ 


We remark that the statement and the proof of Theorem 12.41 do not involve the properties of 
the loss function. In other words, as long as the resulting drive function satisfies conditions (i) and 
(ii), the delay region is bounded. We notice that if the delay region is bounded, then there exists 
a constant s such that {L a > 0} C [0, T] x (0, s). We will utilize Theorem 12.41 repeatedly when we 
discuss the liquidation strategies in subsequent sections. 


2.2 Applications to GBM and Exponential OU Underlyings 

Henceforth, we shall investigate analytically and numerically the optimal liquidation timing when 
S follows (i) the geometric Brownian motion (GBM) model with s) = fi and a(t, s) = a > 0, 
as well as (ii) the exponential OU model with /x(t, s) = f5{6 — log(s)) and a(t, s) = o > 0. 

We will study the liquidation timing of a stock, European put and call options. For both the 
GBM and exponential OU cases, the risk-neutral measure Q is uniquely defined by (12.21) . and the 





Novikov condition is satisfied (see Appendix[A|. Furthermore, the Q dynamics of S is a GBM with 
drift r and the no-arbitrage prices (see 12.311 of a call and a put with strike K and maturity T are 
given by 


C(t, s) = s $(di) - Ke-^ T - t ^{d 2 ), and P(t, s ) = Ke-^ T - t ^(-d 2 ) - s $(-di), (2.14) 

where 3> is the standard normal c.d.f. and 

log (£) + (r + 4 )(T-t) 


di = 


a 


VT^t 


d 2 = d\ — o\fr — t. 


In order to numerically compute the non-trivial liquidation strategy, we solve the variational 
inequality (VI) of the form 


min< — Lf — s)s L“ — 


cr 2 (t, s)s 2 


L“ + rL a - G a , L a } = 0, 


(2.15) 


with terminal condition L a (T,s ) = 0 and where (t,s) € [0, T) x M + . In Secti on [71 we sh ow 
t hat the above VI admits a unique strong solution in the terminology of Bensoussan and Lions 
( 19781 ) under conditions that include the GBM and exponential OU cases (see Theorem 17.21) . For 
implementation, we adopt the Crank-Nicholson scheme for the VI (12.151) on a fi nite (discretized) 
grid D = [smin, S mn ] x [0, T]. We refer to the book by Glowinski ( Glowinski . 19841 . Chap.3) for 
details on numerical methods for solving inhomogeneous Vis of parabolic type. 


3 Optimal Liquidation with a GBM Underlying 

We begin our first series of illustrative examples under the GBM model. In view of Proposition 
we observe that, if p, < r, it is never optimal to hold a stock or a call (or in general any 
positive delta position) regardless whether we introduce a risk penalty or not. On the other hand, 
Proposition 12.21 also implies that, if p > r and a = 0, it is always optimal to delay. 

However, with a non-zero risk penalty (a > 0), the solution can be non-trivial. To see this, we 
note that the drive function associated with a call is given by G^ all (t, s) = {n — r)sC s — mjj((rn — 
C(t, s)) + ) where C(t, s ) is the call price in (I2.14|) . In particular, the penalty term is strictly positive 
at s = 0 and decreasing, and it vanishes for large s. On the other hand, the first term (fi — r)sC s 
is strictly increasing from zero at s = 0. This implies that there exists a price level s such that 
Gc aU (t, s ) is positive in [0, T] x [s,oo). In turn, it follows from (12.121) that the sell region must be 
bounded (possibly empty) and the delay region is unbounded. The same argument applies to the 
case with a stock. Figure [2] illustrates this. 

Next, we consider the liquidation of a put option. Recall the put price P(t,s ) given in (12.141) . 
Its negative Delta implies that for ji > r the drive function Gp ut (t, s ) < 0, V(t, s ), meaning that it 
is optimal to sell immediately by Proposition 12.21 In contrast, when /x < r, the sell region is empty 
if a = 0 , but under risk penalization the optimal strategy may be non-trivial. 


Proposition 3.1. Consider the optimal liquidation of a put under the GBM model with /a < r and 
a > 0. Then, the delay region is hounded. Furthermore, it is non-empty if m < K and 3 1 £ [0, T] 
such that — P(t, 0)) + ) < i — fi. 


9 

























| G>o| | Delay Region | 

- liquidation boundary 

- G = 0 


| G<0 1 | Sell Region | 


Figure 2: The optimal liquidation boundaries (solid) and the zero contours of G a (dashed) of a stock (left 
panel) and a call option (right panel). We take T = 0.5, r = 0.03, /j = 0.08, a = 0.3, K = 50, a = 0.1. The 
loss function is given by %/>(£) = £, with the benchmark m = 50 for the stock and m = (7(0, K) for the call. 


Proof. The drive function for the put G a = Gp ut (t,s ) = (r — fx)s$(—di) — aip((m — P(t,s )) + ) 
satisfies 


lim G a (t, s) = — aijj((m — K) + ) < 0, (3.1) 

s-> 0 

lim G a (t,s ) = < 0, (3.2) 

s—> oo 

dG a 

s) = [r-fi- - P(t, s)) + )l {m>P(M)} ]<h(-(ii) - (r - /j)sT, (3.3) 

where T = > 0. In turn, we fix any b € (0 ,aif(m)) and define if(£) := min{ S}. This 

implies the inequality 


G a (t, s ) := (r — /i)s$(— d\) — — P(t , s)) + ) > G a (t, s). 

Then by Corollary 12.31 we only need to show that G° satisfies the assumptions of Theorem 12.41 
We observe that G is bounded above and it follows from (13.2D that Hindoo G (t, s) — >• —ab < 0 
for every t € [0, T]. Moreover, there exists s > 0 such that for every s > s, — P(t,s)) + ) = b. 

Consequently, we have 


dG° 

~df 


- r)s4>{di) 


log (#)-(r + 4 )(r-t) 
2a(T — t) 2 


< 0 , 


for s > max{s, I\ exp((r + a 2 /2)T)} and t € [0,T], Since G a (0,s) —>• —ab as s —oo, we can 
choose b € (0 ,ab) such that 3k > max{s, K exp((r + <t 2 /2 )T)} and —b > G(0, s) > G(t,s ) in 
[0,T] x [A:,oo). Therefore, G satisfies the assumptions of Theorem 12.41 

Finally, suppose 3 1 € [0, T] such that aif'((m — P(t, 0)) + ) < r — //, where m < K. It follows 
from m and (13.31) that G a (t, 0) = 0 and ^-(t, 0) > 0, so that the set {G a > 0} is non-empty. 
In turn, the inclusion (12. 12|1 implies that the delay region is also non-empty. □ 

Remark 3.2. As an example, the delay region is empty if — P(t, 0)) + ) > r — ^>0, 

[0, T]. Indeed, since we have G a (t, 0) < 0 and 0) < 0, Vi € [0, T\, G a cannot be strictly 

positive. By Proposition \2.Sl it is optimal to sell immediately. 


10 



























Proposition 13.II is illustrated in Figure [3j In these examples, the delay region is non-empty and 
the sell region is unbounded but may be disconnected (Figure [3] (right)). This can arise when, 
for example, G a (t, 0) < 0 for every t € [0, T], but mint max G a (t, s) > 0. The intuition for a 
disconnected sell region is as follows. If the put is deeply in the money (i.e. when St is close 
to zero), its market price has very limited room to increase since it is bounded above Ke~ r ^ T ~ t \ 
At the same time, delaying sale further will incur a penalty. Therefore, when the penalization 
coefficient a is high, it is optimal to sell at a low stock price level. On the other hand, if the put 
is deep out of the money (i.e. when St is very high), the market price and the Delta of the put 
are close to zero, meaning the drive function becomes more negative and selling immediately is 
optimal. 




Figure 3: The optimal liquidation boundary (solid) and the zero contour of G a (dashed) of a put option 
under GBM dynamics with the loss function ip(£) = l. We take m = 2A", a = 0.001 (left panel), and 
m = P(0, AT), a = 0.01 (right panel). Parameters: T = 0.5, r = 0.03, g = 0.02, a = 0.3, K = 50. 


For a long position in calls and the underlying stock, or in puts, the Delta C s takes a constant 
sign. As an example of a derivative with a Delta of non-constant sign, we consider a long straddle. 
This is a combination of a call and a put with strike prices K\ < K 2 respectively and the same 
maturity T. The payoff of a straddle is given by h STD (ST) '■= (St — K\) + + ( K 2 — St) + ■ The 
market price of a long straddle, denoted by C STD , is simply the sum of the respective Black-Scholes 
call and put prices, i.e. C STD (t , s ) = C(t, s ) + P(t, s ). For simplicity, we set K\ = K 2 = K. 


Proposition 3.3. For the optimal liquidation of a long straddle position under the GBM model, it 
follows that 

(i) if i-L — r, the delay region must he empty; 

(ii) if n> r, the delay region is unbounded; 

(in) if n <r, the delay region is bounded. 

Proof. The straddle’s drive function is Gg TD (t, s) = (/x — r)sCg TD (t, s ) — aif((m — C STD (t,s)) + ). 
For p. = r, the conclusion follows immediately by Proposition 12.21 If g, > r we simply notice that 
Gg TD (t, s) —>• 00 as s —> 00 for every t € [0, T], and the assertion follows from the inclusion (12.121) . 

Now suppose [j, < r. We will show that G a satisfies the assumptions of Theorem 12.41 Clearly, 
Gg T D bounded above. Since C STD (t,s ) — >• 00 as s —>• 00 for every t € [0, T], then there exists 
s > 0 such that, for every s > s and t € [0, T], ip((m — C STD (t, s)) + ) = 0. Moreover, for 
s > max{s, K exp ((r + a 2 /2 )T)}, we have 


d$(di) 

dt 


<t>(di) 


log (£)-(r + 4)CT-*) 

2a(T-t)i 


> 0 , 
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and thus 


dG 


STD 


dt 


(t,s) = 2 (fJL-r)s 


d$(di) 

dt 


< 0 . 


This implies Gg TD (0, s) > Gg TD (t,s) for every t € [0, T]. Since Gg TD (0, s) —>• —oo as s —>• oo, for 
a fixed b > 0 there exists kb > 0 such that Gg TD (0,s) < — b for every s > k^. Therefore, setting 
k = max{s, K exp ((r + ct 2 /2)T) , we have Gg TD (t,s) < Gg TD (0, s) < —b in [0,T] x [fc,oo). 
Therefore, the assumptions of Theorem 12.41 are satisfied and we conclude. □ 


In particular, Proposition 13.31 suggests that when p. < r, the sell region is unbounded, even if 
a = 0. In Figured! we illustrate the optimal liquidation boundaries for cases (ii) and (iii). When 
the investor is bullish (left panel: /i = 0.08 > 0.03 = r), the liquidation boundary is increasing 
and the delay region is on top of the sell region. Interestingly, the opposite is observed when the 
investor is bearish (right panel: [i = 0.02 < 0.03 = r). 




Figure 4: Optimal liquidation boundary and the zero contour of G a for a straddle under the GBM model 
with the loss function = i. We set K = 50, m = C STD (0, K), a = 0.1, r = 0.03, ^ = 0.08 (left panel) 
and n = 0.02 (right panel). 

We end this section by discussing the liquidation timing of a stock with an infinite horizon 
(T = oo). This leads to the following stationary optimal stopping problem 

L(s ) = sup E. s | J e~ ru G a (S u )du j . (3.4) 

where G a (s) = (ji — r)s — — s) + ) and T is the set of F-stopping times taking values in [0, oo]. 

When fi < r, selling immediately is optimal according to Proposition 12.21 as for the case with finite 
maturity. As it turns out, the liquidation problem has the opposite trivial solution when /r > r, 
that is, it is optimal to hold forever. 


Proposition 3.4. If ji > r, then the value function L(s ) in (13.41) is infinite and it is optimal to 
never sell the stock. 
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Proof. Consider a candidate stopping time r = oo. Then, by applying Tonelli’s theorem, we have 


C poo h r poo poo 

E s |y e~ ru G a (S u )du'j = E s <j/ e -ru (/i — r)S u du ~ a J e“ r X(m - S u ) + )duj 

poo poo 

— / e (/i-r)u(^ _ — a e~ ru E {^((m — S' u ) + )} du 

JO Jo 

poo pOO 

> / e^~ r J u (^ — r)sdu — a / e~ ru if{m)du = oo, 

JO JO 

since n > r and if} is increasing. Hence, L(s) = oo and it is never optimal to sell. 


□ 


4 Optimal Liquidation with an Exponential OU Underlying 

In the exponential OU model, the stock price satisfies the SDE 

dS t =/3(9-log S t )S t dt + aS t dW t , (4.1) 

with 0 G R and /?, a > 0. Therefore, the optimal liquidation premium L(t, s ) is given by equation 
(12.101) with the drive function 

G a {t,s ) = [/3(9 - log(s)) - r]sV s (t,s ) - - V(t,s)) + ), (4.2) 

where V(t,s ) is a generic option price in (12.31) . 

In contrast to the GBM case, the optimal liquidation strategy can now be non-trivial for a stock 
or a call when there is no penalty. More generally, we can prove that the delay region is in fact 
bounded. The intuition should be clear: when St is very high, it is expected to revert back to its 
long-term mean, so that selling immediately becomes optimal. 

Proposition 4.1. Under the exponential OU model, the delay region for a call is bounded. 

Proof. The drive function G^ all for the call is given by (14.21) with V(t,s) = C(t,s ) (see (|2.14D 
for the call price). It is bounded above, so it satisfies condition (i) of Theorem (12.41) . As is well 
known, the call price satisfies c>C gf s ' > < 0. In addition, /3(9 — log(s)) — r < 0 iff s > exp (9 — -|), and 
> 0 for s > K exp ((r + cr 2 /2)T). In turn, we have 

dG g t al1 (*» s ) = - lo g( s )) - r \ s - G(t , s)) + )l{ m >c(M)I ^ °> 

for s > max{exp(0 — ^),K exp ((r + a 2 /2)T)} and t G [0, T], This implies Gq o11 (0, s ) > Gfg all {t, s). 
Fix b > 0. Since Gf, all (0, s) -G —oo, 3kb > 0 s.t., Vs > kb, G^ all (0,s) < —b. Hence, if we set k = 
max{exp(# — U, K exp ((r + a 2 /2)T) ,kb}, we are guaranteed that Gfj all (t,s) < Gg, aZZ (0,s) < —b 
in [0,T] x [/c,oo), thus satisfying condition (ii) of Theorem 12.41 As a result, Theorem 12.41 applies 
and gives the boundedness of the delay region for a call. □ 

Since a stock can be viewed as a call with strike K = 0, Proposition [6] also applies to the optimal 
liquidation of a stock over a finite time horizon. Also, we notice the delay region can be empty, 
and we can identify this case by finding the maximum of the drive function. As an example, we 


13 











consider the case of the stock with penalty function — St) + ) = (m— St) + , and we obtain the 
maximizer of G a in different scenarios 


arg max G a 


'exp(0 - 1 - 
< exp(6» - 1 - £) 


if exp(0 - 1 - < m, 

if exp(# — 1 — > m, 

otherwise, 


and the corresponding maximum values 


max G a 


((3 — a)§i — a{m — s*) 

< Ph 

k (/3(0 — log(m)) — r)m 


if si < m, 
if §2 > m, 
otherwise, 


where 

r — a v 

Si = exp(0 - 1-—), s 2 = exp(0 - 1 - -). 

Thus, the delay region is non-empty if and only if maxG Q > 0. 

The optimal liquidation boundary for stock is shown in Figure [5] for a = 0 (left panel) and 
a > 0 (right panel). We notice that, in both cases, the optimal strategy is to sell immediately if 
St is high enough. Intuitively, if St is high, it is expected to revert back to its long-term mean, 
so selling immediately becomes optimal. However, if St is low, the optimal behavior depends on 
the parameter a. On one hand, St is expected to increase and thus the investor should wait to 
sell at a better price (Figure [5j right panel). On the other hand, such benefit is countered (if the 
penalization coefficient is high enough) by the risk incurred from holding the position, and this 
induces the investor to sell immediately. As a consequence, the sell region is disconnected (Figure 
[5l right panel). 




Figure 5: The liquidation boundary (solid) and the zero contour of G a (dashed) for a stock under exponential 
OU dynamics. Parameters: T = 0.5, r = 0.03, 9 = log(60), j3 = 4, a = 0.3, = £, a = 0 (left), a = 1.5 

(right). 

Figure [6] illustrates the delay region for a call option with penalty. In the right panel, we 
observe the interesting phenomena where the sell region is connected and contains the nonempty 
delay region. If the parameter /3 (which measures the speed of mean reversion) is not sufficiently 
high, there may be no time for the price of the option to revert back to its long-term mean before 
expiration, so that selling immediately becomes optimal close to maturity. 
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Figure 6: The liquidation boundary (red solid) and the zero contour of G a (dashed) for a call under 
exponential OU dynamics. We take a = 0.2, 6 = log(60), (3 = 4 in the left panel, and a = 0.001, 9 = log(50) 
and (3 = 0.2 in the right panel, with common parameters T = 0.5, r = 0.03, a = 0.3, ip(£) = £. 


Proposition 4.2. For the liquidation of a put option under the exponential OU model, the delay 
region is bounded if and only if a > 0. 

Proof. The drive function is given by 

Gputfa s) = [r- (3(0 - log(s))]s<F(—di) - aip((m - P(t, s)) + ), (4.3) 

If a = 0, then we have {Gp ut > 0} = {s > exp(^ — 0)}. By (12.1211 . the delay region contains this 
set, so it is unbounded. 

Now let a > 0, and we have the limit 


lim Gp ut (t,s) = —aip(m) < 0. 

5—^00 


(4.4) 


Next, we fix any b G (0, aif(m)) and define ip(£) := min b}. With this, we have 
G a (t, s) := [r - (5(9 - log(s))]s4>(—di) - aif {(m - P(t, s)) + ) > Gp ut (t, s). 


We observe that G a is bounded above and by (14.41) lim^oo G a (t, s) —>■ — ab < 0 for every t G [0, T]. 
Moreover, there exists s > 0 such that for every s > s, ip((m — P(t, s)) + ) = b. As a result, we have 


8G a 

~dT 


(/3(9 - log(s)) 


x , x lo s(i 
- r)s(j>(di) -i 1 - 


-(r + 4 )(T-f) 
2 a(T — t) 2 


< 0 , 


for s > max{s, exp(^ — 9),K exp((r + c 2 /2)T)} and t € [0, T], Also, we notice that G a ( 0, s) — > —ab 

as s — >■ —oo. This allows us to choose a b G (0, ab), then there exists k > max{s, K exp((r+a 2 /2)T)} 
such that —b > G(0, s) > G(t , s) in for (t, s ) G [0, T] x [k, oo). Therefore, G satisfies the assumptions 
of Theorem 12.41 By Corollary 12.31 we conclude the boundedness of the delay region. □ 


Proposition 14.21 is illustrated in Figure [71 When St is low, it is expected to revert back to the 
(higher) long-term mean, and the put price will decrease. This generates an incentive to sell at 
a low stock price level. If a = 0, when St is high, there is no reason to sell since the put price 
is very low and expected to increase. Consequently, the delay region is on top of the sell region 
(Figure [71 left panel). However, this is no longer true when we incorporate a non-zero risk penalty 
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| G>0 | Delay Region | 

- liquidation boundary 

- G = 0 


| G<01 | Sell Region ] 


Figure 7: The liquidation boundary (solid) and the zero contour of G a (dashed) for a put option under the 
exponential OU model. We take a = 0 and K = 50 in the left panel, and a = 0.01 and K = 40 in the right 
panel. Common parameters: T = 0.5, r = 0.03, a = 0.3, /3 = 4 and 9 = log(60), 


which reduces the value of waiting. As a result, the holder may sell the put at high and low stock 
prices. In fact, if the penalization coefficient is large and/or when the time-to-maturity is very 
short, the optimal liquidation premium may be zero at all stock price levels, resulting in an empty 
delay region (Figure 0 right panel). 

5 Quadratic Penalty 

As a variation to the shortfall-based penalty, we consider a risk penalty based on the realized 
variance of the option price process from the starting time up to the liquidation time. Precisely, 
the investor now faces the penalized optimal stopping problem 

J a (t, s) := sup E m ( e~ r ^V(r, S T ) -a T e~ r ^d[V, V] u 
r&Tt,r l Jt 

= sup E M ( e~ r ^V(r, S T ) -a [ T e~ r ^a 2 (u, S u )S 2 u V s 2 (u, S u )du} , 
reT t ,T l Jt J 

where [V,V] denotes the quadratic variation of option price process V defined in (12..31 b Figure 
HI illustrates the realized quadratic penalty associated with a simulated call option price path. 
Compared to the shortfall penalty in Figure 12.41 the realized quadratic penalty is increasing at all 
times, even when the option price is above its initial price. 

Following (12.61) . we define the optimal liquidation premium by L a (t,s ) := J a (t,s ) — V(t,s). 
Again, we shall discuss the stock or option liquidation problems under the GBM and exponential 
OU models. 

5.1 Optimal Timing to Sell a Stock 

We first consider the liquidation of a stock with the GBM dynamics in terms of the perpetual 
optimal stopping problem: 


L a (s) 


snpE S y\- ru G a (s u )d U y 


(5.1) 
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Figure 8: Realized quadratic penalty (dashed) based on a simulated price path (solid) of a call under the 
GBM model with a = 0.05. The price path and other parameters are the same as Figure [TJ 


with the drive function G a (s) := (/i — r)s — aa 2 s 2 . If p < r, then selling immediately is always 
optimal since G a is always negative. In contrast if fi > r, then we obtain a non-trivial closed-form 
solution. 


Theorem 5.1. Let /i > r. The value function L a (s ) in ()5.1I) is given by the formula 


L a {s) — { ^ ^ x s x — s + -Bs"! l{ s <s*}, 


2 - A 


where 


B = 


s = 


aa 


2 /r + a 2 — r ’ 
1 — A 


A = 4 

< 7 Z 


a 


o 


— - n + \ — - fi) + 2ra 2 


(2 — X)B ’ 

and the stopping time t* = inf{t > 0 : St > s*} is optimal for dSU). 
Proof. We first show that ()5.21) is the solution of 


(TV 


min <j rA(s) — /isA'(s)-— A"(s) — G a (s), A(s) f = 0, s > 0, 


(5.2) 


(5.3) 

(5.4) 


with A(0) = 0. To do this, we split M + into two regions: V\ = (0, s*) and V 2 = [s*, 00 ) with s* > 0 
to be determined. We conjecture that A(s) = 0 in T> 2 , and for s € T>\ A(s) solves 

2 2 

rA(s) — nsA'(s) — —^—A "(s) — G a (s) = 0. (5.5) 

By direct substitution, the general solution to equation (15.51) is of the form 

A(s) = Cis Al + C 2 S X2 — s + Bs 2 , 
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where C\ and C 2 are constants to be determined, B is specified in (15.31) and 


A k = 


1 


<7- 


<7 


— H + (- 1 ) 


a- 




+ 2 rcr 2 


, fee { 1 , 2 }. 


We apply the continuity and smooth pasting conditions at s = 0 and s = s* to get 


limA(s) = 0 C 1 = 0 , 

40 

lim A(s) = 0 =}► C 2 (s *) A2 - s* + B(s*) 2 = 0, 
si's* 

lim A'(s) = 0 => A 2 C , 2 (s *) A2_1 - 1 + 2 Bs* = 0 . 
sts* 


(5.6) 

(5.7) 


Solving the system of equations (15.61) (15.71) gives C 2 and s* as in (I5.3|) ~ (15.41) . One can verify by 
substitution that A(s) is indeed a classical solution of f[5.5D . 

By Ito’s formula and (15.5p . ( A(St))t>o is a (P,F)-supermartingale, so for every F-stopping time 
r and n € N, we have 

( srAn 'i 


A (s) > E 0 , s |j e~ ru G a (S u )du j . (5.8) 

Maximizing (15.8(1 over r and n yields that A(s) > L"(s) for s > 0. The reverse inequality is deduced 
from the probabilistic representation A(s) = Eo, s | fy e~ ru G a (S u )du ^, with the candidate stopping 
time t* := inf{f > 0 : St > s*}. Hence, we conclude that A(s) = L a (s ) and t* is optimal. □ 


The optimal liquidation threshold s* in (15.41) is non-negative if and only if A 2 < 1, which is 
equivalent to the condition /r > r in Theorem 15.11 Otherwise, L a (s ) = 0 and the optimal strategy 
is to sell immediately. 

In Figure [9} we illustrate the optimal liquidation premium L a (s ) for various values of ^ and a. 
As fjj increases, the optimal threshold as well as the optimal liquidation premium (at all stock price 
levels) increase (left panel). On the other hand, a higher volatility reduces the optimal liquidation 
premium at every initial stock price. We also observe that L a (s ) smooth-pastes the level 0 at the 
optimal threshold s*, as is expected from (15.61) and (15.71) . 

If S follows the exponential OU dynamics, the drive function for liquidating a stock is 

G"(s) = [/3(0 — log(s)) — r — as]s. (5.9) 

In this case, we do not have a closed-form solution. Nevertheless we observe from (15.91) that the 
delay region is non-empty, namely, {L a > 0} A {s < s}, where s is determined uniquely from the 
equation 

/3(0 — log(s)) — r — as = 0. 

On the other hand, since G a —>■ —00 as s —>• 00 , we expect intuitively that the investor will sell 
when the stock price is high. 


5.2 Liquidation of Options 

We now discuss some numerical examples to demonstrate the liquidation strategies for European 
call and put options. With strike K and maturity T, the drive functions are respectively given by 

Gcaiii 1 ^) = s$(di)(^-r--ao- 2 s4>(di)), (5.10) 

Gp ut (t, s ) = s<b(— d\) (r — fi — acr 2 sd>(—di)). (5.11) 
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Figure 9: The optimal liquidation premium for a stock under the GBM model for different values of /j and 
a. In the left panel, we take r = 0.03, a = 0.3 and a = 0.2, and the liquidation threshold s* = 9.37, 7.97, 6.52 
for /i = 0.09,0.08,0.07 respectively. In the right panel, we take r = 0.03, /r = 0.08, and a = 0.1, and the 
liquidation threshold s* = 10.63,7.97,6.26 for a = 0.25,0.30,0.35. 

When /i < r and a > 0, the drive function Gc all (t, s ) is negative for all ( t , s), so it is optimal to 
sell the call immediately. However, when ^ > r and a > 0, we notice from (15.101) that, when the 
stock price is sufficiently large (resp. small), the drive function of a call is negative (resp. positive). 
Hence, as we see in Figure flOl it is optimal to sell the call when the stock price is high, and the 
optimal liquidation boundary is lower as the penalization coefficient increases. In contrast to the 
shortfall penalty, the investor now is subject to a higher penalty when the stock price is high under 
the quadratic penalty. Consequently, the sell region is now above the delay region, as opposed to 
being at the bottom in the shortfall case in Figure [2] (right panel). 

In the put option case, we observe from (|5.11l) that 

lim r — n — aa 2 s<&(—di) = lim r — — ota 2 s&(—d\) = r — /r. 

s—>-0 s—>oo 

Consequently, when /r < r and the stock price is sufficiently large or small, the drive function 
is strictly positive and it is optimal to hold the position. In contrast, the shortfall converges to 
'ijj(m) > 0 as s increases (see (|3.2D ). which means that it is optimal to sell when the stock price 
is high (see Figure [3]). We illustrate the timing strategies under quadratic penalty in Figure [TO] 
(right). As expected there is a low and a high delay regions which are separated by a sell region in 
the middle. Also we notice that as the penalization coefficient a increases, the sell region expands. 

Under the exponential OU model, the drive functions for selling a call and a put are, respectively, 

Gcaiiitis) = S<^{d 1 )(9-r-/3logs-aa^s^idi)), 

Gp ut (t , s) = s<3?(— di) — 9 + (3 log s — aa 2 s&(—di) 

In Figure fill we can visualize the optimal liquidation premium L a {t , s ) for a call (right panel) and 
a put (left panel). In the call case, the delay region, which corresponds to the area where L a > 0, 
is bounded. When s is sufficiently high, L a vanishes and it is optimal to sell. This is intuitive since 
lims^oo Gc al i{t, s) = —oo and G^ all is positive for sufficiently small s. 

In contrast, the drive function for the put Gp ut (t , s ) is negative when s < exp , for every 

a > 0. Therefore, as Figure [IT] indicates, one expects the optimal liquidation premium to vanish for 
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Figure 10: The liquidation boundaries for a call option (left panel) and a put option (right panel) under 
the GBM mode with different values of a. Parameters: T = 0.5, r = 0.03, a = 0.3, K = 50, p = 0.08 (call) 
and \x = 0.02 (put). 


small s, so the investor will sell when the put price is high. Compared to Figure [3 with a shortfall 
penalty, the investor does not sell when the underlying stock price is very high. This is because 
the drive function Gp ut (t, s ) stays positive for large s (recall (I2.12D ). As time approaches maturity, 
the delay liquidation premium decreases to its terminal condition of value zero. 
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Figure 11: The optimal liquidation premium for a call option (left) and a put option (right) with exponential 
OU dynamics. We take T = 0.5, r = 0.03, a = 0.3, K = 50, a = 0.1, /? = 4 and 9 = log(60). 


6 Conluding Remarks 

In summary, we have provided a flexible mathematical model for the optimal liquidation of option 
positions under a path-dependent penalty. We have identified the situations where the optimal 
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timing is trivial, and solved for non-trivial liquidation strategy via variational inequality. The 
penalty type as well as the penalization coefficient can give rise to very different liquidation timing. 
Our findings are useful for both individual and institutional investors who use options for speculative 
investments or risk management purposes. 

For future research, a natural direction is to adapt our model to the problem of sequentially 
buying and selling an option. Moreover, one can cons ider appl ying th e met hodology to derivatives 
other than equity options. For example, we refer to Leung and Liu ( 2012h for a recent study on 
the liquidation of credit derivatives with pricing measure discrepancy but without risk penalty. 
It would be both mathematically interesting and challenging to study option liquidation under 
incomplete markets. On th e other hand, ou r model can be extended to m a rkets with liquidity cost 


and p rice impact (see e.g. lAlmgren ( 200,11 ) : Lorenz and Almgren ( 2011 ): Schied and Schoneborn 
(2003)). Finally, the path-dependent risk penalization can also be incorporated to dynamic portfolio 
optimization problems to account for adverse performance during the investment horizon. 


7 Strong Solution to the Inhomogeneous Variational Inequality 

In this section, we follow the terminology and procedures in Bensoussan and Lionsl ( 1978h . and 
establish the existence and uniqueness of a strong solution to the variational inequality (J2.15D 
under conditions that are applicable to the GBM and exponential OU models. 

Preliminaries. We express prices in logarithmic scale by setting Xt = log (St). Equation (12.11) 
then becomes 

dx t = v (t, x t )dt + »(t, x t )dw t , (7.i) 

for some functions n(t,x) and r](t,x). Next, we define the operator A by 


A[-] 


k 2 (£, x ) d 2 ■ d ■ 


a_ 

dx 


a 2 (t,x 


d_ 

dx 


. . <9 • 

+ a i( t ,x) — +r ’’ 


(7.2) 


where 

ai(t,x) = ■ —/r(/..r) - r](t,x), 

In term of log-prices, we express the drive function as g(t, x) = G a (t, e x ) and the optimal liquidation 
premium as u(t, x) = L(t, e x ). Throughout, we denote the domain V = [0, T] x R. In order to solve 
the VI (12.151) . it is equivalent to solve the VI problem: 


a 2 (t,x ) = 


n 2 (t, x) 


— + A[u] — g(t , x) > 0, u(t, x) > 0, (t, x) € V, 

' {~m + -AH-g(t,x)) u = o, (t,x)£V, (7.3) 

^u(T, x) = 0, x € R. 


We describe an appropriate class of solutions for (17.31) in a suitable Sobolev space and prove 
that such a solution exists and is unique. First, let us define, for A(x) = exp(—n|x|), n € N, 

£|(R) = {u | VXv € L 2 ( R)}, 

Hl(R) = {v€ L 2 m I — € L\( R)}, 

«J,#) = I , lim v{x) = 0}. 

\x\—>oo 
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These are Hilbert spaces when endowed with the following inner products 


(/, 5 )z, 2 = f Xfgdx, f,g£L{(R), 

JR 

( f,g)m = J R x f gdx + ^ e 


We denote by XL], A (K) the set of functions w € "H A (M) with compact support. For u € XLq a 
w € 7^ A (R), we define the operator 

f / f)qj f)qn r) \ \ C f)qj f 

l\(t,u,w) = / a 2 (t, x) ( A— — + w— — ) dx + / ai(t, x)X—wdx + r / Xuwdx. 

J R \ OX OX OX OX J J j[£ OX J 

We can assume without loss of generality (IBensoussan and Lionsl . Il978l . Sect. 3.2.17) that I\ is 
coercive on 77^ A (M), i.e. 


T\(t,w,w) > a:||u;||#i Vw € XL]. A (R), a > 0. 

Integrating by parts allows us to extend T\ to a bilinear form on the whole space Xi] x (M). In 
particular, we set 


X\(t, u, v) 


, ..dudv . .dXdu 
a 2 {t,x) A—— + a 2 (t,x) — —v 
ox ox Ox ox 


dx + 


{ 1 da\ 

V “ 2 ¥ 


A uvdx 


with u,v € nl tX { , 

Following Chapter 5.9.2 of lEvansI (119981 ) and Chapter 2.6 of IBensoussan and Lionsl (119781 ). we 
define the space £ p (0,T;X) consisting of all strongly measurable functions y '■ [0,T] —> X with 


|x||£P(0,T;X) — (J 


\ 1/p 

p x dt\ , 1 < p < oo, 


and for p = oo, 


llxll£°°(o,T;jn = ess sup ||x(t)||x, 

0 <t<T 


For x € T 1 (0, T; X ), we say v € T 1 (0, T; X) is the weak derivative of y, denoted by v = if 

^X{t)dt = -J^ w{t)u(t)dt , VwGC“([0,T]). 

The Sobolev space 'H 1 (0, T;X) consists of all functions y G £ 2 (0,T;X) such that the weak deriva¬ 
tive exists and belongs to £ 2 (0,T;X). Furthermore, we set 


ixlwHo ,T-x) - ^ llxWllx + \\-Q^x(t)\\x dt 


1/2 


(7.4) 


which makes Xi 1 (0,T\X) an Hilbert space (see Chapter 5.9.2 in Evand ( 19981 )). 

Main Results. 
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Definition 7.1. A function u : T> —>■ R is a strong solution of problem (17.311 if,\/v€'H A V > 0 
a.e., the following conditions are satisfied: 

u € T 2 (0,T;^ A (R)), f € £ 2 (0,T;T 2 (R)), 

- (fjM ~ u ) -Z\{t;u,v - u) < (g,v - u), 
u > 0 a.e. in P, 
n(T, x) = 0, x € R. 

We shall impose the following conditions on 02 , 01 , 5 . 

Assumption A. 02 , and ^ € £°°(P); ai and € C°(P); g € 'H 1 (0, T; £|(R)). 

Theorem 7.2. Under Assumption A, the variational inequality in m has a unique strong solu¬ 
tion. 


Proof. Assumption A i s equ ivalent to assumptions (2.223), (2.224), (2.238), (2.239), (2.240) of 
( Bensoussan and Liona . 1978 . Chap. 3), and we also follow their Remark 2.24 to use A(x) = e _n ^ 
for some arbitrarily fixed n > 0 in our definition of Hilbert spaces. In turn, we can apply their 
Theorem 2.21 and our statement follows. □ 


Our main objective is to verify that Assumption A is satisfied for our applications so that 
Theorem 17.21 applies to ensure the existence of a unique strong solution to the VI (|2.15jl . To see 
this, we first write down the operators associated with the log-price X t = log(S)) under the GBM 
and exponential OU models, namely, 


A[v] 


a 2 d 2 v dv 
2 dx 2 ^ dx ’ 


u 2 d 2 v ~ dv 

5? + <»-*>*■ 


Therefore, 02 = cr 2 /2 is constant and a± is an affine function in x for both cases, so these coefficients 
meet the requirements in Assumption A. 

It remains to verify that the drive function g(t,x) = G(t,e x ) € 'H 1 (0, T; £ 2 (M)). In view of 
CZ3D, we want to show that there exists n > 0 such that 



[ [ f<?(C x)e aM'j dxdt, 

Jo Jr ^ ' 


and 



dg_ 

dt 


(t,x) 



2 

dx dt 


are finite, where 

9 (t,x) = (r- n(t,e x ))e x V s (t,e x ) - - V(t,e x )) + ), 

|| (t, x ) = (r - n(t, e x ))e x V ts (t , e x ) + - V ( t , e x )) + )V t (t, e x )l {m>v(t ^ )} . 

Here, the subscripts of V indicate the partial derivatives in t and s. Recall the drift functions 
n(t, e x ) = n under GBM and fi(t, e x ) = (3(9 — x) under exponential OU models. We notice that, in 
both cases, the drift does not depend on t, so we just write ti(e x ). Also, we observe that if and ff are 
increasing, and is bounded for any finite t. For both call and put options, there exist positive 
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constants hi, qi, h 2 , q 2 such that \V s (t,e x )\ < 1, \V t {t, e x )\ < h\e x + qi, | e^)| < h 2 e x + < 72 - 

Together, these imply the time-independent bounds for both models: 

\g(t,x)\ < \r — g(e x )\e x + ai/}(m) = o(e 2 ^), 

\%(t,x)\<\r- fi(e x )\(h ie x + Ql )e x + a^(m)(h 2 e x + q 2 ) = o(e 2 ' x '). 
at 

This implies that by choosing n > 4, we have 
f T f T 

y o lg(t,x)lci(R) dt ^ y o \\\ r ~ v( eX )\ eX + a ^ m )fc\(ji) dt < °°> 

J 0 l||^ a; )ll Cl(R) dt - J o lll r ~ g{e x )\(h 1 e x + qi)e x + a^\m)(h 2 e x + q 2 )\\ 2 c2 ^dt < oo. 

Hence, we conclude that g G 'H 1 (0, T;£|(R)) for both puts and calls under the GBM and exponen¬ 
tial OU models, and Assumption A is satisfied. 


As a final remark, Sect. 3.4 of Bensoussan and Lionsl ( 1978l l also provides the probabilistic 
representation of the strong solution u(t,x ) of the VI (17.31) . given by 


u(t, x ) 


sup E t ,x{ [ e r(u Q g(u, X u ) du\ , 
T&Tt,T Ut J 


(7.6) 


where dX u = rj(u, X u )du + k(u, X u )dW u and X t = x. By the definition L(t,e x ) = u(t,x), the 
optimal stopping problem in (17.61) resembles that for the optimal liquidation premium in (12.101) . 


A Novikov Condition 


Under the GBM model, the Sharpe ratio A = tl-eL is constant, so the Nov i kov condition is clearl y 
met. Let us consider the exponential OU case. By Corollary 3.5.14 of iKaratzas and Shre ve (T991J), 
it suffices to show that for every t € [0,T], 


Be > 0 s.t. E < exp 


\~£ 


2 (u,S u )du ) } < oo. 


(A.l) 


Using Jensen’s inequality and Tonelli’s Theorem, we have, for every e > 0, 


E 


jexpQf \\{u,S u ) 2 dv^ <\ 6 E jei A 6 i ’ S “) 2 j du 

= - E |e2{^ l[e - lo g (s " ) ] 2 + K 2 + ^ 3[ ®- log(s " ) ]} J du. 


for some real constants K\, K 2 and K 3 . Moreover, Cauchy-Schwarz inequality implies that 

E| e f{^ 1 [ 0 - log ( s “)l 2 + /<2 + E ' 3 [ e - log ( s “)]}J < ^/e { e U^i[0-iog(>S tl )] 2 +A' 2 }| E | e£ K 3 (0-iog(5 u ))|. (a.2) 

Since log(S,,) has the normal distribution, there exists e > 0 s.t. E {exp (eK\ [8 — log(S'u)] 2 + K 2 ^} 
and E{exp (eK^O — log(£ u )])} in (IA.2D are finite. Hence, condition (jA.ll) holds. 
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